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Abstract 



We present thermal conductivity measurements on very pure and dense bulk 
samples, as indicated by residual resistivity values as low as 0.5 mW cm 
and thermal conductivity values higher than 200 W/mK. In the normal state 
we found that the Wiedemann Franz law, in its generalized form, works well 
suggesting that phonons do not contribute to the heat transport. The thermal 
conductivity in the superconducting state has been analysed by using a two- 
gap model. Thank to the large gap anisotropy we were able to evaluate 
quantitatively intraband scattering relaxation times of tt and a bands, which 
depend on the disorder in different way; namely, as the disorder increases, it 
reduces more effectively the relaxation times of tt than of a bands, as suggested 
by a recent calculation . 

l.INTRODUCTION 

Since the discovery of MgB 2 [0, experimental findings established a phonon mediated 
s-wave superconductivity. Despite its standard origin, superconductivity in MgB 2 showed 
several unusual properties that can be ascribed to the presence of two gaps with different 
amplitudes emphasized by tunnelling and specific heat measurements. Theoretical studies 
@ pointed out that the peculiar electronic structure is the origin of this behaviour, 
being the larger gap, A a , associated with two-dimensional a bands and the smaller one, A n , 
associated with three-dimensional 7r bands. The two-gap model offers a simple explanation 
of several anomalies in the superconducting state || 0. Multiband effects due to the 
different parity of a and ir bands were predicted in the normal state, too [0, and, recently, 
have found confirmation in transport measurements ||. Since interband impurity scattering 
turns out to be negligible, different bands behave as separate conduction channels in parallel 
and either a or ir channel prevails, depending on the disorder degree, being the clean (dirty) 
samples dominated by 7r(<r) conduction. 

In this paper we perform a quantitative analysis of thermal conductivity measurements 
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in the superconducting state to investigate the role of a and it bands depending on the 
sample purity. 

Thermal conductivity k in the superconducting state among the other transport prop- 
erties gives information on quasi-particle (QP) excitations and their dynamics with the 
advantage of probing only the QP response, since the superffuid does not carry heat. On 
the other hand, a major complication in the analysis of the thermal conductivity is often 
a substantial phonon contribution to the heat current. Consequently, the interpretation of 
experimental data can be ambiguous. 

From a basic point of view, QP condensation below Tc causes a decreasing of the electron 
contribution to the thermal conductivity, K e , and an increasing of the phonon contribution, 
n p ; thus k below Tc can show a shoulder or a peak depending whether the heat current is 
dominated by electrons or phonons, and a more complex behaviour is exhibited when both 
the contribution are important. 

From the beginning the thermal conductivity of MgB 2 showed the unexpected feature 
of not exhibiting any signature of the superconducting transition [[J [TO] [[□]] ]T2|. This 



"anomaly" was initially ascribed to a perfect compensation of n e and n p [TO] or to the 



presence of large thermal resistance at the grain boundaries [IT]. As a matter of fact, the 
discovery that MgB 2 presents two gaps, one of which is so small that QP condensation 
becomes exponential only below a reduced temperature of the order of t = T /Tc ~ 0.2, 
opened new perspectives in the interpretation of thermal conductivity data. Hence, the two 
gap model successfully used to fit specific heat data can also be applied to the thermal 
conductivity to estimate the fraction of energy carried by the two bands, providing a useful 
mean to investigate multiband effects. 

2. EXPERIMENTAL DETAILS 

Dense (up to 2.4 g/cm 3 , 90% of the theoretical density), clean and hard cylinder shaped 
samples have been prepared by a single step method [TB| similar to the one reported in ref 



|T4| |15| . Amorphous or crystalline B and Mg, put in Ta crucibles welded under argon and 
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closed in quartz tubes under vacuum, were heated up to 950°C . An X rays spectrum of a 
sample prepared by crystalline B is shown in fig. 1. All the peaks associated with the MgB 2 
phase are present, no extra peaks due to the presence of free Mg , MgO, are detected. In 
the inset, SEM (Scanning Electron Microscopy) image of the same sample is shown. The 
image shows a network of well connected grains (2-4 mm large). 

Two specimens were prepared for physical measurements: one from crystalline boron 
(MGB-1S) , one by using enriched 11 B (MGB11-1S). The samples were cut in the shape of 
parallelepiped bar (1x2-3x12 mm 3 ). 

The thermal conductivity was measured using a steady state flux method with a heat 
flux sinusoidally modulated at low frequency (u=0. 003-0. 01 Hz ). Under these conditions the 
thermal conductivity is extracted as k = J(z/)/VT(z/), where J(u) is the heat flow provided 
at the frequency v and VT(z/) is the temperature gradient oscillating at the frequency v. 
A small resistive heater (lxl mm 2 ) is glued by GE varnish on the top end of the bar, 
being the bottom of the bar thermally connected with the sample holder. In such way a 
longitudinal heat flow is assured. The gradient applied to the sample was varied from 0.1 
to 0.3 K/cm. Seebeck effect was measured simultaneously with the thermal conductivity 
providing a precise determination of the critical temperature. 

3. RESULTS AND DISCUSSION 

Resistivity measurements of the samples from 30 to 300 K are shown in fig. 2. In table 
1 we report T™ set , the amplitude of the transition ATc, p(40K) and the residual resistivity 
ratio defined as RRR = p(300K) / p(40K) for the two samples. The excellent quality of the 
samples is proved by the high T c values (38.9, 38.7 K), the small AT C values (0.2, 0.3 K), 
the low values of p(40K) (0.6, 2.5 pQcm) and the large values of RRR (7-15). Between the 
two samples the enriched 11 B (MGB11-1S) has the lower resistivity values which can mainly 
be ascribed to the very good quality of the Eagle-Picher enriched n B 

Thermal conductivity measurements of the two samples from 4 K to 250 K are shown 
in fig. 3. The outstanding quality and high density of the samples is evident from thermal 
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transport properties as well. Indeed, these samples have more than one order of magni- 



tude higher thermal conductivities than p oh/crystalline samples || [|Kj ||TT|. In particular 
MGB11-1S exhibits a thermal conductivity as large as 215 W/Km at 65 K which is nearly 
two times higher than those of a single crystal |12| , proving the excellent purity and density 
of this sample. 

The thermal conductivity of MGB11-1S increases monotonically in the superconducting 
state, a change of slope is observable at about 8 K, while no signature of the superconducting 
transition is present; the normal state curve exhibits a pronounced maximum at about 70 
K. Similar behaviour is presented by MGB-1S even if only data above 10 K are available. 

The normal state behaviour of both samples is typical of a good metal in which the 
electron contribution to heat transport prevails. In order to estimate the relative weight 
of K e and K p we can consider the effective Lorenz number defined as L e ff = np/T. This 
quantity is assumed equal to Lq = 2.45 x 10~ 8 WQK~ 2 by the Wiedemann- Franz law (WFL); 
in good metals L e ff is equal to L at low temperature where the electron impurity scattering 
prevails, then it decreases showing a minimum at about one tenth the Debye temperature 
G>£), deeper and deeper with increasing sample purity JIB]] ; in dilute alloys where k p is not 



at all negligible L e ff becomes larger than L and the ratio L e ff/L ~ (1 + K v j ' n e ) gives 
information on the relative weight of K e and n p ||17|| . In fig. 4 we plot L e ff for the two 



samples from 40 to 250 K. The curves exhibit the typical behaviour of good metals with 
slightly different levels of purity: at low temperatures the curves approach L from below; 
they show a minimum at around 130 K (~ O.IB^, considering that in MgB 2 Qd ~1000 
K [|18j [pl|]) which is more pronounced for MGB11-1S; finally, the curves increase towards 
Lq. Thus, we can conclude that in the normal state the WFL substantially works and 
hip can be neglected. Similar results were obtained also in sintered sample with low thermal 
conductivity ||, while different conclusions have been drawn on a single crystal ||12|| , where 



a violation of the WFL was claimed. Actually, in a small sample with not well defined 
geometrical shape as a single crystal the geometrical factor which relates conductivity to 
conductance can be different in thermal and electrical measurement, giving an uncorrect 



evaluation of the Lorenz number. In large bar cut specimens the geometrical factors can be 
estimated with better precision and more reliable verification of the WFL can be made. 

In the superconducting state, since k p decreases as (T/0d) 3 , we can assume that the 
thermal conductivity is dominated by electrons as well. Only at very low temperature 
T <^Tc due to QP condensation which decreases n e and enhances k p , their relative weight 
can change. Hence, we assume that lattice vibrations give negligible contribution to the 
thermal transport in the temperature region of our interest and in the following we analyse 
the thermal conductivity data in term of the electron contribution only. 

The electron thermal conductivity in the superconducting state can be written as: 

K s e {T) = K n e {T)g{t,a) (1) 

where k™(T) is the electron thermal conductivity in the normal state and g(t,a), for 
a given reduced gap, a = A(0) / KTc, and reduced temperature, t = T/Tc , takes into 
account the QP condensation. The function g(t,a) was calculated in the framework of 
the BCS theory in the dirty and clean limit cases [[2(J |21[] reproducing very well thermal 
conductivity in the low temperature superconductors. To analyse experimental data using 
eq. (1) two main problems have to be solved. First, eq. (1) has to be generalized to 
multiband conduction; second, k™ (T) has to be estimated. We start from the latter issue. 

At low temperature, if the scattering with impurities prevails, k™(T) can be obtained by 
the WFL, k" (T) = LqT/ po where po is the residual resitivity. But in clean samples scattering 
with phonons has also to be considered. In ref. || this was done using a generalized WFL. 
We can write: 

W e = W\ + Wl (2) 

where, for the Matthiessen's rule, the thermal resistance W e = l//c e is the sum of the 
thermal resistivity for scattering with impurities, Wl, and for scattering with phonons, Wf. 



Following ref. |T6[ we can write: 



w: = -p- (3) 
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P P (T) 
WF = 



2/3 



L T 

,5 „«, 



(4) 



where p p (T) = Ap'Q D Iq D ^ T ^ e x_i^^L e -x^ is given by the Bloch-Gruneisen equation, 

pi is the temperature coefficient and n a is the average number of electrons per unit cell in a 
given band. 

Despite its simplicity, this model describes very well the normal state thermal conduc- 
tivity of samples with very different degrees of disorder, with reasonable parameter values 
||; it contains the right features, vanishing to zero at low temperature linearly with T and 
requiring only three free parameters (p , 6^ and C = ^p'® D -^ (jf) ^ ) whose reliability 
can be checked by resistivity . The best fitting procedure is performed from 40 K to 200 
K and the curves obtained with the parameter values listed in tab. 2 are reported in fig. 4 
as continuous lines. The theoretical curves fit the experimental data in the normal state in 
excellent way, while just below Tc, they lay above the data, decreasing linearly with tem- 
perature, where the experimental data start to decrease with larger slopes, due to the QP 
condensation. In tab. 2 we can see that the Gd and C values, which determine the intrinsic 
term , are nearly the same for both the samples, while the po values change by a factor 4, 
as consequence of the different purity of the two samples; moreover we can see that the po 
values are slightly lower than p(40K) reported in table 1. Anyway, po, Op and C values are 
in fair agreement with those obtained by fitting resistivity measurements ||. 

Now we address the problem of generalizing eq. (1) in the case of multiband conduction. 
Following the approach used for the specific heat ||, we can write: 

-f-^r = xg(t, a n ) + (1 - x)g(t, a a ) (5) 

K e K 1 ) 

where a v = A„-(0) / KT C and o a = A (T (0) / KT C ; the relative weights x and (1 —x) that in 
ref. [H are related to the energy which condenses in each band, in our case rather represent 
the energy fractions carried by the it and a bands, respectively; thus the relative weights 
take into account the mobility of the carriers in each band, also. 
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If we indicate with k™(T) = xk" (T) and k™(T) = (1 — x)k™ (T) the thermal conductivity 
in the normal state of 7r and a bands, eq. (5) states that the tt and a bands conduction 
occurs in parallel, as suggested by ref . [|TJ . 

Once a suitable g function is considered eq. (5) allows to fit the thermal conductivity 
data in the superconducting state with x, a n and a a as free parameters. To describe a and 



tt QP we choose the g function in the dirty limit |20j in agreement with the fact that below 
40 K fv"(T) is dominated by the scattering with impurities. 

In figure 5 we show k/T in the superconducting state for the two samples. The contin- 
uous line represents «™(T)/T, the dashed line represents K s e {T)/T given by eq. (5). The 
parameters obtained with the best fit procedure are summarized in table 2. We can see that 
the two-gap model provides an excellent agreement with the experimental data. We find 
^=0.57 ( 0.6) and a CT =2.17 (1.9) for MGB11-1S (MGB-1S); these values, which correspond 
to A(0) 7r = 1.9(2.0) meV and A(0) 7r = 7.2(6.3) meV are in agreement with those found 
by specific heat data cr n ~0.6-0.65 and o a ~1. 9-2.2. For the x parameter we find 0.85 and 
0.75 for MGB11-1S and MGB-1S, respectively. From the specific heat analysis it comes out 
that the energy fraction which condenses in each band is roughly the same; thus our results 
imply that carriers in tt bands are more mobile than in a bands. The lacking signature 
of superconducting transition in thermal conductivity, naturally follows from the two-gap 
model; in fact, we find that the main contribution to heat transport comes from carriers in tt 
bands whose strong condensation starts for cr n /t » 1 which means t < 0.2. In MGB11-1S 
data a change of slope occurs at about 8 K which in this framework represents the tt QP 
condensation; in fact in figure 5 the dotted line represents K s e (T)/T given by eq. (5) with 
x = 1 and a n =0.57 . The curve fits well the data just below 10 K, while from 10 to 40 K 
also the a QP excitations contribute to the transport. 

The right weight of tt and a contributions to transport can be calculated for both samples 
by k%(T) = XKg(T) and k£(T) = (1 — x)k™(T) once the reliability of the x and Kg(T) 
evaluation has been verified. We point out that in the fit procedure, varying slightly the gap 
values, the quality of the fit does not change very much, while the x parameter for both the 

8 



samples is very well defined. On the other hand, below Tc, k"(T) = L^T/po and therefore 
it only depends on the residual resistivity p which can be experimentally estimated. So we 
can write: 

<(T) = L T— = LoT<£ < (6) 
Po 

<(T) = L T— = L^le.rl (7) 
Po 

where is the dielectric constant, uj v ^ a are the plasma frequencies and r^ CT are the 
intraband scattering relaxation times with impurities for tt and a bands. 

Our fits state that carriers in tt bands mainly contribute to carry heat. Indeed, the 
plasma frequency is larger for tt than for a bands || and this is enhanced in polycrystalline 
samples by averaging in the three direction. On the other hand, the scattering relaxation 
times change from sample to sample. Introducing in eq.s (6) and (7) the values of x and po 
listed in tab. 2 and uj P7V = 6.226 eV and co> P(T =3.403 eV || we obtain the scattering relaxation 
times for the two samples summarized in tab. 3. For MGB11-1S we find -7^=2.2 xlO -13 s 
and r*=1.3xl0 -13 s and the residual mean free paths l na .=r^ a v F7TtCT , where vp^^ are the 
averaged Fermi velocities for it an a bands (wi? 7r =5.6xl0 5 m/s and Vp a =3.2xl0 5 m/s ||), 
come out 1^=1.2 xl0 _7 m and / CT =4.1 xl0~ 8 rn. These values are very large, exceeding the 
lattice constants by more than two orders of magnitude, indicating the excellent purity of 
this sample. For MGB-1S the relaxation rates are reduced, but the residual mean free paths 
are still large (7 7r =2.9 xlCT 8 m and l a =1.8 xlO -8 m). 

Now we can look at the relaxation times in more details. In MGB11-1S t\ >t^, while 
in MGB-1S both the relaxation times are lower and <r l a . Practically, as the disorder 
increases, it reduces more effectively r\ than r x a . Actually, it is not easy to say which kind 
of disorder is present in MGB-1S, being itself a quite pure sample. But typical defects in 
MgB 2 are vacancies or substitutions in the Mg site, which form more easily than in the B 
site. In these cases the relaxation rate for intraband impurity scattering is larger in tt than 
in a bands (r^ <T l a ) |T|. This prediction, which was well verified for large amount of defects 
||, is here confirmed also going from very pure to rather pure samples. 
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This result can open new perspectives in thermal transport properties. In fact, as the 
disorder increases, we expect that the transport of a bands prevails, and in superconducting 
state, the rapid QP condensation due to the large gap related to a bands can become evident. 
Thus, in disordered samples, we expect that the thermal conductivity would diminish in 
absolute values, but it should show a wide shoulder below T c . This fact, if true, is quite 
unusual, in fact the disorder generally smoothens, rather than enhancing features. 

As a matter of fact, looking at the thermal conductivity data in literature, polycrystalline 
samples do not show more evident shoulder than clean samples. We think that polycrys- 
talline samples, that surely are more disordered than single crystals or bulk samples, present 
resistive grain boundaries, which contribute to the thermal resistance masking the intrinsic 



behaviour of superconducting grains [TT| . To emphasize the transport of carriers in a bands 
by the progressive inhibition of transport in 7r bands, it is necessary to gradually introduce 
defects in pure samples. This can be done by suitable chemical substitution with the advan- 
tage of introducing a controlled amount of defects in a chosen site. Up to now, substituted 
samples have been obtained by sintering powders previously synthesized from the pure el- 
ements and such samples are not suitable for transport measurements. A second way is to 
introduce disorder by irradiation, but in this case the problem is to obtain a uniform defect 
distribution. 

In conclusion, the role of disorder in MgB 2 , which multiband effects make so peculiar, 
has been studied in pure samples by thermal transport measurements. We analysed in detail 
the thermal conductivity in the superconducting state, and thank to the large gap anisotropy 
we were able to evaluate quantitatively the intraband scattering relaxation times of 7r and 
a bands. 
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Figure captions 

Figure 1. X-Rays pattern diffraction of an enriched 11 B bulk sample. The inset shows a 
SEM image of the same sample. The length scale of the picture is indicated in the bottom. 
Figure 2. Resistivity measurements from 40 to 300 K. 

Figure 3. Thermal conductivity measurements of the samples from 10 K to 250 K: 
the best fitting curves obtained with the parameter values listed in tab. 2 are reported as 
continuous lines. 

Figure 4. The effective Lorenz number L e ff/Lo = np/(TLo) from 40 to 250 K. 

Figure 5. k/T in the superconducting state. The continuous line represents the calcu- 
lated k™/T] the dashed line represents n s e /T given by eq. (5) with the parameter values 
summarized in tab. 2; the dotted line represents K s e /T given by eq. (5) with x = 1 and a n 
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given in tab. 2. 



Table captions 

Table 1. Critical temperature T ( 2/ lse ', amplitude of the transition ATc, p(A0K) and 
residual resistivity ratio defined as RRR = p(300K) / p(40K) . 

Table 2. Values of the parameters 6, po and C obtained from the normal state thermal 
conductivity best fit and a n , a a and x obtained by fitting eq. (5) with the superconducting 
state thermal conductivity. 

Table 3. The intraband scattering time with impurities for tc and a bands, r£ and r^and 
the residual mean free paths l n}a =r^. a Vp^^ a (v Fn =5.6xl0 5 m/s and v Fa =?).2xlQ b m/s are 
the averaged Fermi velocities for 7r an a bands ||). 



Table 1 



sample 


T° nset , K 


AT C ,K 


p(A0K),pQcm 


RRR 


MGB11-1S 


38.7 


0.2 


0.58 


15.3 


MGB-1S 


38.9 


0.3 


2.1 


7.1 



Table 2 



sample 


9,K 


Po, p^lcm 


C.pVtcmK 2 






X 


MGB11-1S 


1190 


0.50 


12 


0.57 


2.17 


0.85 


MGB-1S 


1130 


1.9 


18 


0.60 


1.90 


0.75 



Table 3 



sample 




4,s 


l n , m 


l a ,m 


MGB11-1S 


2.2xl0- 13 


1.3xl0 -13 


1.2xl0 -7 


4.1xl0~ 8 


MGB-1S 


0.5xl0~ 13 


0.6xl0~ 13 


2.9xl0~ 8 


1.8xl0 -8 
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2e [degrees] 




T(K) 



